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measuring at least one effect of the change on the incident ultrasound; 
using the measured at least one effect to determine an absorption coeffi- 
cient for the radiation in the region; and using the determined absorption 
coefficient to determine concentration of the component in the region. 
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PHOTOACOUSTIC ASSAY AND IMAGING SYSTEM 

FIELD OF THE INVENTION 

The invention relates to non-invasive in-vivo methods and apparatus for determining 
the concentration of a substance in a body and for determining the concentration of the 
5 substance as a function of position in the body. 

BACKGROUND OF THE INVENTION 
Methods and apparatus for in-vivo and in-vitro measurements of blood glucose levels 
are known in the art. Generally, the methods and apparatus are relatively complicated and 
measurements of a person's blood glucose levels are usually performed in a clinic or 

1 0 laboratory with the aid of a technician and costs of the measurements are relatively high. 

Methods and apparatus for determining blood glucose levels for use in the home, for 
example by a diabetic who must monitor blood glucose levels frequently, are available. These 
methods and associated devices are generally invasive and usually involve taking blood 
samples by finger pricking. Finger pricking is perceived as inconvenient and unpleasant and 

15 to avoid finger pricking diabetics tend to monitor their glucose levels less frequently than is 
advisable. Moreover, many conventional glucometers require routine purchasing of sample 
sticks and pricking needles, which is bothersome and adds cost to the user. There is a need for 
glucometers that are easy to use and that perform non-invasive in-vivo assays of blood sugar. 

PCT Publication WO 98/38904, the disclosure of which is incorporated herein by 

20 reference, describes a "non-invasive, in-vivo glucometer" that uses a photoacoustic effect in 
which light energy is converted to acoustic energy to measure a person's blood glucose. Pulses 
of light at a wavelength for which light is absorbed by glucose is directed by the glucometer 
to illuminate a part of the person's body, such as a fingertip, comprising soft tissue. The light 
pulses are typically focused to a relatively small focal region inside the body part and light 

25 from the light pulses is absorbed by glucose and converted to kinetic energy in a region of 
tissue in the neighborhood of the focal region. The kinetic energy causes temperature and 
pressure of the absorbing tissue region to increase and generates acoustic waves, hereinafter 
referred to as "photoacoustic waves", that radiate out from the absorbing tissue region. An 
acoustic sensor comprised in the glucometer contacts the body part and senses the 

30 photoacoustic waves. Intensity of the waves is a function of the concentration of glucose in 
the absorbing tissue region and their intensity as measured by the sensor is used to assay the 
glucose. 
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However, light is scattered by body tissue and even though the light is focused to a 
small focal region inside the body, the location and size of the absorbing tissue region are not 
accurately known. As a result, the generated photoacoustic effect and measurements of the 
person's glucose levels are not necessarily the result only of glucose concentration in the 
5 person's blood. Characteristics of the absorbing tissue region, such as density of blood vessels 
therein, that can affect concentration of glucose in the absorbing region are often not 
accurately known. Measurements of blood glucose levels can therefore be affected by 
unknown variables that substantially compromise the reliability of the measurements. 

. US Patent 5,941,821 describes another non-invasive in-vivo glucometer that uses a 

10 photoacoustic effect to assay blood glucose. Light at a wavelength at which glucose absorbs 
light is modulated at a suitable frequency and directed by the glucometer to illuminate a 
region of a person's body. Glucose in blood and interstitial fluid in tissue near the surface of 
the region absorbs the light and converts the absorbed energy to kinetic energy that heats the 
tissue. Temperature of the tissue increases and decreases cyclically in cadence with the 

15 modulation of the light. The alternate heating and cooling of the tissue results in periodic 
heating of air in contact with the surface of the illuminated region, which generates sound 
waves in the air. A microphone comprised in the glucometer provides measurements of 
intensity of the sound waves that are used to determine a concentration of glucose. 

SUMMARY OF THE INVENTION 

20 An aspect of some embodiments of the present invention relates to providing a non- 

invasive, in-vivo glucometer that determines a person's glucose level substantially only from 
glucose concentration in the person's blood. 

An aspect of some embodiments of the present invention relates to providing a 
glucometer that locates at least one blood vessel in a person's body using ultrasound. The 

25 glucometer determines the person's glucose level by assaying glucose present substantially 
only in blood in the located at least one blood vessel. 

In embodiments of the present invention, the glucometer comprises at least one 
ultrasound transducer that radiates ultrasound into the person's body and at least one 
ultrasound sensor that receives energy from the radiated ultrasound that is reflected or 

30 transmitted by features in the body. The glucometer locates the at least one blood vessel from 
the reflected/transmitted energy using methods known in the art. In some embodiments of the 
present invention, the glucometer comprises at least one light source that illuminates a region 
of tissue in which the at least one blood vessel is located with at least one pulse of light at a 
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wavelength for which light is absorbed by glucose. Glucose in a volume of blood, hereinafter 
referred to as a "bolus" of blood, in the at least one blood vessel absorbs energy from the at 
least one light pulse. The absorbed energy generates changes in the bolus that are a function of 
the amount of energy absorbed, which in turn is a function, inter alia, of the concentration of 
5 glucose in the bolus. 

At least one change in an acoustic property of the bolus and/or acoustic phenomena 
generated in the bolus by the absorbed energy is measured by the glucometer and used to 
assay glucose in the bolus. In some embodiments of the present invention a measurement of a 
change and/or phenomenon is performed during illumination of the bolus with the at least one 
10 light pulse. In some embodiments of the present invention a measurement is performed after 
illumination of the bolus with the at least one light pulse. 

According to an aspect of some embodiments of the present invention a photoacoustic 
wave generated by the absorbed energy is used to assay glucose in the bolus. 

When the region of tissue in which the at least one blood vessel is located is 
15 illuminated with the at least one light pulse, the bolus and generally tissue surrounding the 
bolus generate photoacoustic waves responsive to energy that they absorb from the at least 
one light pulse. The glucometer senses the photoacoustic waves and identifies locations of 
their origins using methods known in the art. The glucometer compares the identified 
locations with the location of the bolus to identify which of the photoacoustic waves 
20 originates in the bolus. The amplitude of the photoacoustic wave so identified is a function of 
the concentration of glucose in the bolus and is used to assay the person's glucose. 

According to an aspect of some embodiments of the present invention, the glucometer 
measures concentration of glucose in the bolus by reflecting ultrasound waves from the bolus. 

Energy absorbed by the blood bolus from the at least one light pulse generates changes 
25 of temperature and pressure in the bolus that change the acoustic impedance of the bolus. The 
change in the acoustic impedance changes the acoustic reflectivity of the bolus. The 
glucometer transmits ultrasound waves into the person's body that are incident on the bolus. 
Reflections of acoustic energy from the incident ultrasound waves by the bolus axe sensed to 
determine the change in. reflectivity. The measured change in reflectivity is used to assay 
30 glucose in the bolus. 

According to an aspect of some embodiments of the present invention the glucometer 
reflects ultrasound waves from the bolus at times during which a photoacoustic wave is being 
generated by the bolus to measure a change in the bolus from which to determine glucose 
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concentration. During generation of the photoacoustic wave, the surface of the bolus is 
moving and ultrasound waves reflected from the bolus are therefore Doppler shifted. 
Generally, during generation of the photoacoustic wave, surface regions of the bolus move at 
speeds that are substantial fractions of the speed of sound in the illuminated tissue. As a result, 
5 frequencies of ultrasound waves that are reflected by the bolus from ultrasound waves 
incident on the bolus are Doppler shifted substantially compared to the frequency of the 
incident ultrasound waves. Intensities of these reflected waves and magnitude of their Doppler 
shifts are sensed and used to determine concentration of glucose in the bolus. According to an 
aspect of some embodiments of the present invention, the glucometer transmits ultrasound 

10 waves through the bolus to measure concentration of glucose in the bolus. 

Temperature and pressure changes in the bolus that result from absorption of energy 
from the at least one light pulse change the speed of sound in the bolus. The change in the 
speed of sound changes a transit time of ultrasound transmitted through the bolus. The transit 
time of at least one pulse of ultrasound through the bolus is measured before and after 

15 absorption of energy from the at least one light pulse to measure the change of the speed of 
sound in the bolus and thereby the concentration of glucose in the bolus. 

In some embodiments of the present invention transit time of ultrasound through the 
bolus is determined from a pulse of ultrasound energy that repeatedly "bounces" back and 
forth between walls of a blood vessel in which the bolus is located before the energy exits the 

20 blood vessel and the bolus. As a result, the ultrasound energy traverses the bolus a plurality of 
times and speed of sound in the bolus is determined from a transit time of the acoustic energy 
pulse through the bolus over a path length many times a characteristic dimension of the bolus. 
Accuracy of a determination of the speed of sound generally increases as the path length over 
which the transit time is measured increases. A speed of sound in blood determined from a 

25 pulse of ultrasonic energy that traverses the bolus a plurality of times is therefore relatively 
more accurate than a speed of sound determined from a pulse of ultrasound energy that 
traverses the bolus once. 

In should be noted that passing an ultrasound pulse through a blood bolus a plurality 
of times, in accordance with an embodiment of the present invention, can also be used to 

30 provide accurate measurements of other changes in the blood bolus caused by absorption of 
energy from the at least one light pulse. For example, "multipass measurements" can be used 
to acquire accurate measurements of changes in the absorption coefficient in blood for 
ultrasound. 
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According to an aspect of some embodiments of the present invention, the glucometer 
measures a change in the bolus by determining wavelengths of sound for which transmission 
or reflection of ultrasound from the bolus is a maximum or a minimum. 

The transmission of sound through a layer of material is a function of the wavelength 
5 of the sound, thickness of the layer, acoustic impedance of the material in the layer and 
acoustic impedances of material on either side of the layer. For a given thickness of the layer 
and values of the impedances, sound transmission through the layer and reflection from the 
layer have maxima at particular "resonant" wavelengths. 

In accordance with some embodiments of the present invention a resonant 

10 transmission and/or reflection wavelength is determined for a bolus of blood prior to 
illumination with the at least one pulse of light. Following illumination of the bolus with the 
at least one light pulse, the acoustic impedances of the bolus and the tissue on either side of 
the bolus are changed. The change in the impedances changes the resonant wavelength. The 
new resonant wavelength is determined, in accordance with an embodiment of the present 

15 invention, by measuring transmission of ultrasound through, and/or reflection of ultrasound 
from, the bolus at a plurality of wavelengths. The amount by which the new resonant 
wavelength is shifted from the old resonant wavelength is a function of glucose concentration 
in the blood bolus. The magnitude of the shift is used to determine glucose concentration in 
the bolus. It should be noted that the "resonances" become sharper and a magnitude of a shift 

20 in a resonant wavelength becomes more pronounced as the impedance differences between the 
bolus and tissue on either side of the bolus increases. 

In accordance with some embodiments of the present invention a Doppler shift of 
ultrasound transmitted through the bolus is measured to assay glucose in the bolus. 

In such embodiments the glucometer transmits ultrasound through the bolus while the 

25 bolus is exposed to the at least one light pulse. During exposure, as glucose in the bolus 
absorbs energy from the at least one light pulse, temperature of the bolus rises and speed of 
sound in the bolus changes. As a result, frequency of ultrasound exiting the bolus is Doppler 
shifted with respect to the frequency of ultrasound incident on the bolus. The magnitude of the 
Doppler shift can be shown to be equal to f(D/C 2 )(dC/dt) where f is the frequency of the 

30 incident ultrasound, D is the path length of the ultrasound through the bolus, C is the speed of 
sound in the bolus and t is time. The expression for the Doppler shift assumes that the speed 
of sound is constant along the path length through the bolus. The derivative, dC/dt, is 
proportional to the time rate of absorption of energy from the at least one light pulse. The time 
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rate of energy absorption is proportional to the concentration of glucose in the bolus and 
intensity of the at least one light pulse. 

The glucometer has been described as illuminating the bolus with light at a 
wavelength for which light is absorbed by glucose. However, it is not possible to choose a 
5 wavelength for the light for which components of the bolus other than glucose do not absorb 
the light. As a result, an acoustic effect measured by the glucometer resulting from absorption 
of light by the bolus at any wavelength is due to absorption of the light by substances other 
than glucose in the bolus, such as cholesterol, albumin and various lipids, as well as by 
glucose. In order to isolate a contribution to the effect due to glucose, and thereby the 

10 concentration of glucose in the bolus, the bolus is preferably illuminated with pulses of light 
at a plurality of different wavelengths and the effect measured at the different wavelengths. 
Using known absorption cross-sections for light by substances in the blood, the contribution 
to the effect from glucose and thereby the concentration of glucose in the bolus is determined. 
It is readily concluded from remarks made in the previous paragraph that the invention 

15 is not limited to assaying blood glucose. In the process of determining glucose concentration, 
in accordance with an embodiment of the present invention, absorption coefficients of other 
substances in the blood are determined and the absorption coefficients may be used to 
determine concentrations of these substances. For example if "N" different wavelengths of 
light are used to determine glucose concentration, in accordance with an embodiment of the 

20 present invention, absorption coefficients of N different components of the blood bolus are 
determined and may be used to estimate concentrations of the components in the bolus. 
Different particular substances in the bolus whose concentrations are to be estimated can be 
assayed, in accordance with an embodiment of the present invention, by proper tc tuning" of 
frequencies used to illuminate the blood bolus. Furthermore, the present invention is not 

25 limited to assaying blood components. Some embodiments of the present invention can be 
used to assay components of other tissues and features of the body, such as for example, 
components of interstitial fluid, blood clots or plaque in blood vessels that are located using 
ultrasound. 

Whereas only light has been described as being used to introduce energy to a blood 
30 bolus, other forms of electromagnetic radiation, such as microwave or RF radiation, may be 
used to impart energy to components being assayed in a tissue volume, in accordance with an 
embodiment of the present invention. Some embodiments of the present invention are used to 
provide a spatial map of concentration of a substance in a region of the body. In some 
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embodiments, a sub-region of the region of the body for which concentration of the substance 
is to be mapped is illuminated with collimated light that is absorbed by the substance. As 
noted above, because of scattering of light in body tissue, the size and location of the sub- 
region is not accurately known. However, changes in acoustic properties of tissue in the sub- 
5 region can be determined for highly localized 'Voxels" in the sub-region. The changes can be 
determined by sensing photoacoustic waves generated in the voxels and/or characteristics of 
ultrasonic waves that are, reflected from and/or transmitted through the voxels, in accordance 
with embodiments of the present invention. Any of the changes in acoustic properties that 
affect transmitted, reflected or generated acoustic waves, such as those discussed above, can 

10 be localized to voxels in the sub-region having dimensions in a range from a few to about ten 
wavelengths of ultrasound that is used to probe the sub-region. Since these changes are a 
function of concentration of the substance, the concentration of the substance in different 
voxels of the sub-region can be determined, in accordance with preferred embodiments of the 
present invention. A map of the concentration of the substance in the body region is acquired 

15 by moving the illuminated sub-region to different "scan" positions in the body region. 

Since different organs and/or features in the body region will usually be distinguished 
by different levels of concentration of a substance, a concentration map of the substance, 
provided in accordance with an embodiment of the present invention, will, in general, image 
the organs and/or features. As a result some embodiments of the present invention are used to 

20 image organs and or features of organs in a region of a body. For example, LDL cholesterol 
and oxidized LDL cholesterol are highly concentrated in plaque. By mapping concentration of 
LDL cholesterol or oxidized LDL cholesterol in a region of a person's body, images of plaque 
deposits on walls of blood vessels in the region can be acquired. 

There is therefore provided, in accordance with an embodiment of the present 

25 invention, A method for assaying a component of a localized region of interest in a body 
comprising: a) illuminating the region with at least one pulse of radiation having a wavelength 
at which the radiation is absorbed by the component to generate a change in an acoustic 
property of the region; b) transmitting ultrasound so that it is incident on the region; c) 
measuring at least one effect of the change on the incident ultrasound; d) using the measured 

30 at least one effect to determine an absorption coefficient for the radiation in the region; and e) 
using the determined absorption coefficient to determine concentration of the component in 
the region. 
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Optionally illuminating the region with at least one pulse of radiation at at least one 
other wavelength, repeating b-d and using the determined absorption coefficient for the at 
least one other wavelength to determine concentration of the component in the region. 

In some embodiments of the present invention, an effect of the at least one effect 
5 comprises a change in acoustic energy reflectance of the region. 

Optionally measuring the change in reflectance comprises: reflecting at least one first 
pulse of ultrasound from the region before illumination of the region with the at least one 
radiation pulse; detennining a first acoustic energy reflectance from the region using the 
amplitude of the at least one reflected pulse; reflecting at least one second pulse of ultrasound 
10 from the region after illumination of the region with the at least one radiation pulse; 
determining a second acoustic energy reflectance from the region using the amplitude of the 
second at least one reflected pulse; detennining a difference between the first and second 
acoustic energy reflectances. 

Optionally measuring the change in reflectance comprises: measuring acoustic energy 
15 reflectance from the region before illumination with the at least one radiation pulse; 
measuring amplitude of Doppler shifted ultrasound waves reflected from the incident 
ultrasound during illumination of the region with the at least one light pulse; using the 
amplitude of the Doppler shifted ultrasound waves to determine acoustic energy reflectance 
that characterizes the region during illumination with the at least one radiation pulse; and 
20 determining a difference between the reflectance before illumination and during illumination. 

In some embodiments of the present invention, an effect of the at least one effect 
comprises a change in speed of sound in the region. 

Optionally measuring the change in speed of sound comprises: transmitting an 
ultrasound wave from a first transducer towards a second transducer so that the ultrasound 
25 wave traverses the region; illuminating the region while the ultrasound wave is present in the 
region; determining a first transit time from the first to the second transducer of a first portion 
of the ultrasound wave, which first portion is present at a first location in the region while the 
region is illuminated; detennining a second transit time from the first to the second transducer 
of a second portion of the ultrasound wave which second portion is present at a second 
30 location in the region at a time at which the first portion is located at the first location; and 
using the difference in the transit times to determine a change in the speed of sound through 
the region at a location in the region between the first and second locations. 
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In some embodiments of the present invention the region is sandwiched between a 
first and a second acoustic interface at which ultrasound is partially reflected and measuring 
the change in the speed of sound comprises: measuring before illumination a difference 
between a transit time of ultrasound between a first and second location that is reflected from 
5 the first interface and a transit time of ultrasound energy between the first and second 
locations that is reflected from the second interface; determining a speed of sound before 
illumination for the region from the transit time difference; measuring after illumination a 
difference between a transit time of ultrasound between a first and second location that is 
reflected from the first interface and a transit time of ultrasound energy between the first and 

10 second locations that is reflected from the second interface; determining a speed of sound after 
iUumination for the region from the transit time difference; and detennining the change in the 
speed of sound from the determined speeds of sound before and after illumination. 

In some embodiments of the present invention, the region is sandwiched between two 
acoustic interfaces at which ultrasound is partially reflected and measuring the change in the 

15 speed of sound comprises: a) introducing a pulse of ultrasound into the region before 
illumination so that the pulse bounces back and forth between the interfaces and repeatedly 
traverses the region; b) sensing energy from the pulse that exits the region through one of the 
interfaces for each of a plurality of back and forth traversals of the pulse through the region; c) 
determining a speed of sound in the region from the times at which the energy is sensed; d) 

20 introducing a pulse of ultrasound into the region after illumination so that the pulse bounces 
back and forth between the interfaces and repeatedly traverses the region and repeating steps b 
and c; and e) determining a difference between the determined speeds of sound. 

In some embodiments of the present invention, the region is sandwiched between two 
acoustic interfaces at which ultrasound is partially reflected and measuring the change in 

25 speed of sound comprises: transmitting ultrasound which is incident on the region prior to 
illumination of the region with the at least one radiation pulse; determining a first frequency 
for which transmission of the ultrasound through the region is a maximum or a minimum; 
transmitting ultrasound which is incident on the region after illumination of the region with 
the at least one radiation pulse; determining a second frequency for which transmission of 

30 ultrasound through the region is a maximum or a minimum and for which the wavelength of 
propagation of ultrasound in the region for the first and second frequencies is substantially the 
same; using a difference between the first and second frequencies to determine the change in 
the speed of sound in the region. 
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In some embodiments of the present invention, the region is sandwiched between two 
acoustic interfaces at which ultrasound is partially reflected and measuring the change in 
speed of sound comprises: transmitting ultrasound which is incident on the region prior to 
illumination of the region with the at least one radiation pulse; determining a first frequency 
5 for which reflection of the ultrasound through the region is a maximum or a minimum; 
transmitting ultrasound which is incident on the region after illumination of the region with 
the at least one radiation pulse; determining a second frequency for which reflection of 
ultrasound through the region is a maximum or a minimum and for which the wavelength of 
propagation of ultrasound in the region for the first and second frequencies is substantially the 

10 same; using a difference between the first and second frequencies to determine the change in 
the speed of sound in the region. 

In some embodiments of the present invention, an effect of the at least one effect 
comprises a change in the frequency of ultrasound that traverses the region. 

Optionally measuring the change in frequency comprises: transmitting an ultrasound 

15 wave from a first transducer towards a second transducer so that the ultrasound wave traverses 
the region; illuminating the region with the at least one light pulse while the ultrasound wave 
is present in the region and wherein the at least one light pulse has a pulse width substantially 
longer than the transit time of ultrasound through the region; determining a first frequency 
shift in the frequency of a first portion of the ultrasound wave, which first portion is present at 

20 a first location in the region while the region is illuminated; determining a second frequency 
shift in the frequency of a second portion of the ultrasound wave which second portion is 
present at a second location in the region at a time at which the first portion is located at the 
first location; and using a difference in the first and second frequency shifts to determine a 
frequency shift that occurs in the ultrasound at a location in the region between the first and 

25 second locations. 

In some embodiments of the present invention, an effect of the at least one effect 
comprises a change in absorption of ultrasound along a path length in the region. 

Optionally the region is sandwiched between two acoustic interfaces at which 
ultrasound is partially reflected and wherein measuring the change, in absorption over the path 

30 length comprises: a) introducing a pulse of ultrasound into the region before illumination so 
that the pulse bounces back and forth between the interfaces and repeatedly traverses the 
region; b) sensing an amount of energy from the pulse of ultrasound that exits the region 
through one of the interfaces for each of a plurality of times that the pulse travels back and 
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forth through the region; c) determining amounts of energy absorbed from the pulse for each 
back and forth traversal through the region; d) introducing a pulse of ultrasound into the 
region after illumination so that the pulse bounces back and forth between the interfaces and 
repeatedly traverses the region and repeating b and c; and e) detemrining the change in 
5 absorption over the path length using the determined amounts of absorbed energy. 

In some embodiments of the present invention, the method comprises determining a 
location of the region. Optionally detennining a location for the region comprises transmitting 
ultrasound into the body. 

In some embodiments of the present invention, transmitting ultrasound incident on the 
10 region and measuring an effect of the change on the incident ultrasound is mediated by at least 
one ultrasound transducer coupled to the body and wherein the method comprises measuring 
acoustic properties of the coupling between the at least one transducer and the body. 

In some embodiments of the present invention, at least one pulse of radiation 
comprises a microwave pulse. In some embodiments of the present invention, at least one 
15 pulse of radiation comprises a pulse of RP energy. In some embodiments of the present 
invention, the at least one pulse of radiation comprises a pulse of light. 

There is further provided, in accordance with an embodiment of the present invention, 
a method for assaying a component of a localized region of interest in a body comprising: a) 
illuminating the region with at least one pulse of light having a wavelength at which light is 
20 absorbed by the component; b) sensing photoacoustic waves generated responsive to the light 
pulse and determining locations of their origins; c) using intensity of those photoacoustic 
waves having an origin in the region to determine an absorption coefficient for the light in the 
region; d) using the determined absorption coefficient to determine concentration of the 
component in the region. 

25 In some embodiments of the present invention, the at least one light pulse comprises at 

least one train of light pulses radiated at a pulse repetition rate. 

Optionally sensing and determining locations for the photoacoustic waves comprises: 
focussing an acoustic reference beam on the region, which reference beam has a frequency 
that is shifted from the pulse repetition frequency of the at least one light pulse train by an 
30 offset frequency and an intensity that causes the region to respond non-linearly to acoustic 
stimuli; and detecting acoustic waves at the offset frequency. 

In some embodiments of the present invention the method comprises illuminating the 
region with at least one pulse of light at at least one other wavelength, repeating b and c and 
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using the determined absorption coefficient for the at least one other wavelength to determine 
concentration of the component in the region. 

Optionally, for a method in accordance with an embodiment of the present invention, 
for which sensing and determining locations for the photoacoustic waves comprises focussing 
an acoustic reference beam, the method comprises illuminating the region with at least one 
other light pulse train characterized by a different wavelength and pulse repetition rate, 
repeating b and c, and using the determined absorption coefficient for the at least one other 
light pulse train to determine concentration of the component in the region. Optionally, the 
light pulse train and at least one other light pulse train illuminate the region substantially 
simultaneously. 

In some embodiments of the present invention, sensing photoacoustic waves is 
mediated by at least one ultrasound transducer coupled to the body and wherein the method 
comprises measuring acoustic properties of the coupling between the at least one transducer 
and the body. 

In some embodiments of the present invention, the method includes determining a 
location of the region. Optionally determining a location for the region comprises transmitting 
ultrasound into the body. 

Optionally, determining a location for the region comprises locating the region to 
within an uncertainty of less than 5 wavelengths of the ultrasound in at least one direction. 
Optionally, determining a location for the region comprises locating the region to within an 
uncertainty of less than 3 wavelengths of the ultrasound in at least one direction. Optionally, 
determining a location for the region comprises locating the region to within an uncertainty of 
about a wavelength of the ultrasound in at least one direction. 

Optionally, the at least one light pulse comprises a pulse of IR light Optionally, the at 
least one light pulse comprises a pulse of visible light 

In some embodiments of the present invention, the method comprises measuring 
acoustic properties of material in the body. Optionally, a determined absorption coefficient is 
adjusted responsive to the measured acoustic properties. 

In some embodiments of the present invention, the body is a living body. Optionally, 
the living body is a human body. 

In some embodiments of the present invention, the region is a bolus of blood in a 
blood vessel of the body. Optionally, the component is glucose. 
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There is further provided, in accordance with an embodiment of the present invention, 
a method of imaging internal features of a body comprising: assaying a component of the 
features for a plurality of regions in the body according to a method in accordance with an 
embodiment of the present invention; and displaying a result of measurements performed in 
5 assaying the component as a function of location of the regions to provide an image of the 
features. Optionally the method comprises overlaying the image provided from assaying the 
component on another image of the region acquired by a different imaging modality. 
Optionally the different imaging modality comprises ultrasound imaging. In some 
embodiments of the present invention a same at least one ultrasound transmitter or ultrasound 
10 sensor is used to acquire both the image provided from assaying the component and the 
ultrasound image. 

In some embodiments of the present invention, the body is a living body. In some 
embodiments of the present invention, the living body is a human body. 

In some embodiments of the present invention, the feature is a plaque deposit in a 
15 blood vessel. In some embodiments of the present invention, the component is LDL 
cholesterol or oxidized LDL cholesterol. In some embodiments of the present invention, the 
feature is a tumor. 

There is further provided, in accordance with an embodiment of the present invention 
a method for detecting photoacoustic waves generated in a region of a body by light 
20 comprising: illuminating the region with pulses of the light at a pulse repetition frequency to 
generate photoacoustic waves in the region; heterodyning the generated photoacoustic waves 
with an acoustic reference beam having a frequency that is shifted from the pulse repetition 
frequency by an offset frequency; and detecting acoustic waves at the offset frequency. 

Optionally, heterodyning comprises focusing the acoustic reference beam in the region 
25 so that acoustic energy in the region from the reference beam is sufficient to cause the region 
to respond non-linearly to acoustic stimuli. 

BRIEF DESCRIPTION OF FIGURES 

Non-limiting examples of embodiments of the present invention are described below 
30 with reference to figures attached hereto.- In the figures, identical structures, elements or parts 
that appear in more than one figure are generally labeled with the same numeral in all the 
figures in which they appear. Dimensions of components and features shown in the figures are 
chosen for convenience and clarity of presentation and are not necessarily shown to scale. The 
figures are listed below. 
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Figs. 1A-1C schematically show a glucometer measuring a person's blood glucose by 
sensing photoacoustic waves generated in blood in a region of the person's body by 
iUumination of the region with light absorbed by glucose, in accordance with an embodiment 
of the present invention; 

5 Figs. 2A-2C schematically show a glucometer measuring blood glucose by sensing 

ultrasound reflections from blood, in accordance with an embodiment of the present 
invention; 

Figs. 3A and 3B schematically show a glucometer measuring blood glucose by 
sensing changes in velocity of sound in blood, in accordance with an embodiment of the 
10 present invention; 

Figs. 4 schematically shows a glucometer measuring blood glucose by sensing 
Doppler shifts in ultrasonic waves reflected from blood while the blood is being rapidly 
heated by illumination with light absorbed by glucose, in accordance with an embodiment of 
the present invention; and 
15 Figs. 5A-5B schematically show a glucometer measuring blood glucose by sensing 

Doppler shifts in ultrasonic waves transmitted through blood while the blood is being heated 
by illumination with light absorbed by glucose, in accordance with an embodiment of the 
present invention; and 

Figs. 6A-6B schematically show a glucometer assaying and imaging cholesterol in a 
20 person's body, in accordance with an embodiment of the present invention. 

DETAILED DESCRIPTION OF EMBODIMENTS 
Figs. 1A-1C schematically show a glucometer 20 measuring a person's blood glucose 
by sensing photoacoustic waves generated in blood in a part 24 of the person's body, in 
accordance with an embodiment of the present invention. Glucometer 20 comprises at least 
25 one ultrasound transducer 26, a light source 28 that provides light at a wavelength at which 
light is absorbed by glucose and a controller 30 that controls the at least one ultrasound 
transducer and the light source. 

At least one ultrasound transducer 26 is acoustically coupled to region 24 and is used 
to profile acoustic impedance and velocity of sound in region 24 and identify and locate 
30 features in the region using ultrasound imaging techniques and methods known in the art. In 
particular, at least one ultrasound transducer 24 is used to identify and locate at least one 
blood vessel 22 in body part 24 using methods known in the art. For example, location of 
blood vessel 22 can be performed by detecting a Doppler shift in ultrasound reflected from the 
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blood vessel caused by velocity of blood in the blood vessel, assuming that a component of 
the velocity of the blood is parallel to ultrasound incident on the blood vessel. Or, location of 
a blood vessel can be performed by sensing reflections of ultrasound from walls of the blood 
vessel. 

5 Numerous and varied ultrasound transducers, configurations of transducers and 

ultrasound imaging methods suitable for use in practicing the present invention and locating 
blood vessels are known in the art. 

< 

At least one transducer 26 may comprise a single transducer or a plurality of 
transducers. At least one transducer 26 may comprise a relatively small transducer that 

10 performs as a point source when generating ultrasound waves and performs as a point sensor 
that senses ultrasound waves in a relatively large solid angle. At least one transducer 26 may 
comprise a shaped transducer that generates a collimated beam of ultrasound and senses 
ultrasound in a relatively small solid angle. At least one transducer may comprise a phased 
array of ultrasound transducers that can be controlled by controller 30 to generate and steer a 

15 . beam of ultrasound. 

At least one transducer 26 may comprise a piezoelectric transducer or other transducer 
known in the art, such as an "optical" ultrasound transducer that converts pulses of optical 
energy into acoustic energy. Optical ultrasound transducers are described in an article entitled 
"Theory of Detection of Shear Stress Pulses with Laser Picosecond Acoustics", by O. 

20 Matsuda, 11th ICPPP, Japan 2000, the disclosure of which is incorporated herein by 

« * 

reference. Such transducers can provide and detect ultrasound pulses having sub-picosecond 
pulse lengths. 

In some embodiments of the present invention same transducers are used to both 
generate and sense ultrasound waves. In some embodiments of the present invention at least 

25 one transducer that is used to transmit ultrasound is not used to sense ultrasound. 

In Figs 1A - 1C and figures that follow, for convenience of presentation, at least one 
ultrasound transducer 26 is assumed, by way of example, to both transmit and sense 
ultrasound and to be controllable to generate a collimated and/or focused pulse of ultrasound. 
In Fig. 1 A controller 30 controls at least one ultrasound transducer 26 to transmit at least pne 

30 collimated pulse of ultrasound waves into body part 24, which ultrasound pulse is represented 
by arc line segments 32. Some of the energy in at least one ultrasound pulse 32 is reflected 
from a localized region 36 of at least one blood vessel 22 and returns to at least one transducer 
26 as a reflected at least one ultrasound pulse 38. 

*♦ 

4 
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A reflected pulse 38 is sensed by at least one ultrasound transducer 26, which 
generates signals responsive thereto that the ultrasound transducer transmits to controller 30. 
Controller 30 processes the received signals, using methods known in the art, to identify and 
locate at least one blood vessel 22. For example, assume a bolus of blood 40 moving through 
5 region 36 of at least one blood vessel 22 has a component of velocity parallel to pulse 32 or 
reflected pulse 38. Reflected pulse 38 may then be identified as originating in at least one 
blood vessel 22 by sensing a Doppler shift in the reflected pulse resulting from the pulse being 
reflected from bolus 40. The location of blood vessel 22 may be determined from a time lapse 
between a time at which the pulses of ultrasound waves 32 are transmitted and a time at which 

10 Doppler shifted ultrasound pulse 38 arrives at ultrasound transducer 26. In addition, reflected 
pulses 38 may arise from reflections of energy from transmitted pulses 32 by walls of at least 
one blood vessel 22 as a result of the walls having acoustic impedance different from the 
acoustic impedance of tissue adjacent to the walls. 

It should be noted that energy from at least one ultrasound pulse 32 is not only 

15 reflected by at least one blood vessel 22 and features thereof. In general, energy from at least 
one ultrasound pulse 32 is reflected by tissue interfaces and regions in body part 24 in which 
the acoustic impedance is changing rapidly. In some embodiments of the present invention, 
energy in the reflections that is sensed by at least one transducer 26 is used to determine 
acoustic properties, such as for example acoustic impedance and the speed of sound, of 

20 regions in body part 24. In some embodiments of the present invention acoustic energy 
reflections are used to determine acoustic properties of coupling of at least one transducer 26 
to body part 24. Measurements of the acoustic properties of body part 24 and of the coupling 
of at least one transducer 26 to the body part are hereinafter referred to as "acoustic calibration 
measurements". 

25 In some embodiments of the present invention, ultrasound generated by at least one 

transducer 26 has a frequency greater than 5MHz. In some embodiments of the present 
invention the frequency is substantially equal to or greater than 10MHz. Sound propagates in 
body tissue at a velocity of about 1.5mm per microsecond and has a useful penetration depth 
for detection of features in tissue, which is limited by absorption and decreases as frequency 

30 increases. For an ultrasound frequency of about 10MHz, the wavelength of ultrasound 
generated by the transducer is about 0.15 mm and a useful penetration depth of about 35 mm. 
For this frequency, glucometer 20 can determine location of at least one blood vessel 22 in a 
direction along which the transducer transmits ultrasound pulses with an accuracy of about 
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0.15mm (i.e. about a wavelength). In some embodiments of the present invention, transmitted 
ultrasound pulses 32 are collimated so that they have a cross-sectional diameter of about 
1.5mm. Glucometer 20 therefore has spatial resolution and penetration depth sufficient to 
locate blood vessels having dimensions typical of dimensions of blood vessels located, by 
5 way of example, in a person's wrist, forearm or groin. 

In Fig. IB, after at least one blood vessel 22 has been located, controller 30 controls 
light source 28 to illuminate body part 24 with at least one pulse 42 of collimated light 
"aimed" at region 36 of at least one blood vessel 22. As light pulse 42 enters body part 24, 
tissue in body part 24 scatters the light and the collimated light spreads out laterally. 

10 Therefore as light pulse 42 penetrates deeper into body part 24 it illuminates an increasing 
volume of tissue in the body part. In the region of bolus 40 a volume of tissue greater than the 
volume of bolus 40 is illuminated. It is therefore seen that, because of scattering, collimated 
light cannot, generally, be used to illuminate substantially only a relatively small, accurately 
defined region of tissue in body part 24. As a result, photoacoustic waves can be generated by 

15 light pulse 42 in locations in body part 24 that are not restricted to blood bolus 40 and the 
light pulse may generate photoacoustic waves at locations substantially removed from the 
location of the blood bolus. In Fig. IB, lines 44 schematically represent an envelope that 
defines the spatial extent of at least one light pulse 42 perpendicular to the direction of travel 
of the light pulse and spreading of lines 44 inside body part 24 indicate lateral spreading of at 

20 least one light pulse 42. 

It should be noted that whereas light source 28 is shown, by way of example, located 
at a side of ultrasound transducer 26 other positions of light source 28 are possible and can be 
advantageous. In some embodiments of the present invention, for example, ultrasound 
transducer 26 is formed with a hole and light source 28 is positioned to transmit light pulses 

25 through the hole. Transmission of light pulses through a hole in ultrasound transducer 26 
provides relatively uniform and more symmetric illumination of a "field of view" of the 
transducer in body part 24 in which the transducer detects features of the body part. 
Furthermore whereas light source 28 is shown as a single light source, light source 28 may 
comprise a plurality of light sources that provide light at wavelengths suitable for determining 

30 glucose concentration in bolus 40. 

Fig. 1C schematically shows photoacoustic waves, which are represented by sets of 
concentric circles, that are generated in body part 24 as a result of illumination by light 42. An 
innermost concentric circle in a set of concentric circles schematically represents a location of 
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an origin of a photoacoustic wave. Sets of concentric circles 48 drawn with dashed lines 
schematically represent photoacoustic waves that are generated by absorption of energy from 
light pulse 42 at locations outside of blood vessel 22 and bolus 40. Set 50 of concentric circles 
drawn in solid lines schematically represents a photoacoustic wave having an origin in bolus 
5 40. 

Ultrasound transducer 26 senses photoacoustic waves 48 and 50 and transmits signals 
responsive thereto to controller 30. In accordance with an embodiment of the present 
invention controller 30 determines a location for the origin of each of photoacoustic waves 48 
and 50 using results of acoustic calibration measurements performed previously and noted in 

10 the discussion of Fig. 1A. Controller 30 then compares the locations to the location of bolus 
40, which was determined previously, as indicated in Fig. 1A. Methods for determining the 
locations of origins of photoacoustic waves are described in PCT Publication WO 98/14118 
and US Patent 5,713,356, the disclosures of which are incorporated herein be reference. 

Controller 30, in accordance with an embodiment of the present invention, is therefore 

15 able to discriminate between photoacoustic waves that are generated in bolus 40 and 
photoacoustic waves that are generated outside of bolus 40. Controller 30 determines that 
photoacoustic wave 50 is located in bolus 40 and that photoacoustic waves 48 are located 
outside of bolus 40. In accordance with an embodiment of the present invention, controller 30 
therefore uses the locations and amplitudes of photoacoustic waves 50 and 48 and calibration 

20 measurements to determine a concentration of glucose in the person ! s blood. 

The amplitude of pressure from photoacoustic wave 50 and the shape of photoacoustic 
wave 50 that is sensed by ultrasound sensor 26 is a function of an amount of energy absorbed 
by bolus 40 from at least one light pulse 42. The amount of energy absorbed from at least one 
light pulse 42 is, of course, a function of an absorption coefficient in blood for light in light 

25 pulse 42. In general the amount of energy absorbed is relatively small, as a result of which, 
the amount of energy is approximately proportional to the absorption coefficient of light 42 in 
blood. The relationship between the amplitude of a photoacoustic wave and an amount of 
energy absorbed by a region of tissue that generates the photoacoustic wave is described in 
US Patent 4,385,634 to Bowen, the disclosure of which is incorporated herein by reference 

30 and in PCT publication WO 98/14118 referenced above. Expressions for the amplitude of a 
photoacoustic wave are also given in an article by Lai, H. M. and Young, K.J. in Acoust. Soc. 
Am. Vol 76, pg 2000 (1982), in an article by MacKenzie et al., "Advances in Photoacoustic 
Noninvasive Glucose Testing", Clin. Chem. Vol 45, pp 1587-1595 (1999) and in an article by 
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C.G.A. Hoelen et al., "A New Theoretical Approach To Photoacoustic Signal Generation", 
Acoust. Soc. Am. 106 2 (1999) the disclosures of which are incorporated herein hy reference. 

If P represents the amplitude of photoacoustic wave 50 then, adopting the expression 
for amplitude given by Hoelen et al., P = A(pC 2 /H)(I 0 Te" ac l)ab, where A is a constant of 
5 proportionality, P is the temperature expansion coefficient for blood, C the speed of sound and 
H the specific heat capacity of blood. I Q and % are respectively the intensity and pulse length 
of at least one light pulse 42. a is the absorption coefficient for light in at least one light pulse 
42 in tissue of body part 24 and d is a path length that light pulse 42 travels in the body part to 
reach blood bolus 40. is the absorption coefficient of blood in bolus 40, which is a function 

10 of glucose concentration in the blood, for light in at least one light pulse 42. The exponential 
factor e" a d is an attenuation factor by which intensity of at least one light pulse 42 is 
attenuated along the path length d. The attenuation factor may be estimated from d, which is 
known from the location of blood bolus 40 and experimentally known values for a. In some 
embodiments of the present invention, the absorption coefficient a is determined as a function 

15 of distance along the path length d using the equation for P given above or in one of the other 
referenced documents and a finite element analysis of the measured photoacoustic signal. 

The wavelength of light in at least one light pulse 42 is chosen so that the light is 
absorbed by glucose. The absorption coefficient otb is therefore a function of a known 
absorption cross section of glucose for light at the wavelength of light in at least one light 

20 pulse 42 and an unknown concentration of glucose in the blood of bolus 40. A measure of the 
amplitude P of photoacoustic waves generated in blood bolus 40 can therefore be used to 
determine concentration of glucose in blood bolus 40. In some embodiments of the present 
invention, controller 30 adjusts amplitude P using results of the acoustic calibration 
measurements to remove biases in P introduced by acoustic properties of tissue in body part 

25 24 and coupling of at least one transducer 26 to the body part. Controller 30 uses the adjusted 
P to determine glucose concentration. 

Glucose has absorption peaks at a number of different wavelengths that can be used in 
the practice of the present invention. For example, glucose has absorption peaks in the mid 
infrared (IR) at 9.7 microns in the combination region at 2. 10, 2.27, and 2.32 microns, in the 

30 first overtone region at 1.73, 1.69, and 1.61 microns, and in the near infrared, with relatively 
low absorption, in bands centered at 0.76, 0.92* and 1.00 micron that can be used in the 
practice of the present invention. However, there is no wavelength for the light at which only 
glucose in the blood absorbs the light. Many different substances in the blood such as 

19 



t 



i 



WO 02/15776 PCT/IL01/00740 

cholesterol, albumin and various fats and proteins absorb light at or near wavelengths of light 
at which glucose absorbs light As a result, the absorption coefficient is a function not only 
of glucose concentration in blood bolus 40 but also of the concentrations (and absorption 
coefficients) the other absorbing substances in the blood bolus. Therefore, in accordance with 
5 an embodiment of the present invention, to assay the person's blood glucose, absorption 
coefficients of light by a bolus 40 of the person's blood are determined from measured 
amplitudes of photoacoustic waves for a plurality of different wavelengths of light. The 
concentration of glucose in the person's blood is then determined from analysis of the 

■ 

determined absorption coefficients using algorithms known in the art. 

> 

1 0 Appropriate criteria for the choice different wavelengths of light used to determine 

glucose concentration and the methods and techniques for analyzing the absorption 
coefficients are well known. US patent 5,957,841 to Maruo et al., US Patent 5,452,716 to V. 
Clift, US Pat 5,348,002 to Caro, and U.S. Patent 4,975,581 to Robinson et al. the disclosures 
of which are incorporated herein by reference, describe methods of determining glucose 

15 concentration from absorption measurements at a plurality of wavelengths. US Patent 
5,957,841 describes determining glucose concentration in tissue from measurements of 
absorption of light in the tissue in three wavelength bands: 1.48 microns - 1.55 microns, 1.55 
microns - 1.65 microns and 1.65 microns - 1.88 microns. Suitable light sources, such as lasers, 
laser diodes (LDs), and light emitting diodes (LEDs) for providing light at wavelengths used 

20 to assay glucose, in accordance with embodiments of the present invention, are readily 
available. For methods of assaying glucose, in accordance with embodiments of the present 
invention, such as a method discussed below with reference to Fig. 4, in which continuous 
illumination of a blood bolus at desired frequencies is required, filtered lamps can also be 
used for providing light. 

25 In some embodiments of the present invention, photoacoustic waves generated in 

blood bolus 40 by at least one light pulse 42 are detected by heterodyning the photoacoustic 
waves with a reference acoustic wave, in accordance with an embodiment of the present 
invention. 

To illustrate the heterodyning detection technique, in accordance with an embodiment 
30 of the present invention, assume that at least one light pulse 42 comprises a train of light 
pulses radiated at a first frequency. Photoacoustic waves will then be generated by blood 
bolus 40 at substantially the first frequency. Assume that the reference source is controlled to 
generate "reference" acoustic waves at a second frequency, offset from the first frequency, 

20 



WO 02/15776 PCT/IL01/00740 

which are focused on blood bolus 40. Assume further that the reference acoustic waves are 
sufficiently intense so that in the presence of the reference acoustic waves the blood bolus has 
a substantial non-linear response to acoustic stimuli. As a result, photoacoustic waves 
generated in blood bolus 40 by light pulses 42 are a function not only of the intensity and 
5 frequency of the light pulses but are a function also of the intensity and frequency of the 
reference acoustic waves. The "acoustic non-linearity" of blood bolus 40 caused by the 
reference waves couples the reference waves and the photoacoustic response of the blood 
bolus to light pulses 42. 

In particular, if the phase of the reference waves and the envelope of optical energy of 

10 light pulses 42 at blood bolus 40 are stable, the coupling of the reference waves and the 
photoacoustic response of the blood bolus will generate a "heterodyned photoacoustic signal". 
The heterodyned photoacoustic signal will have a frequency substantially equal to the beat 
frequency of the first and second frequencies will be transmitted from the blood bolus. The 
pressure amplitude of the heterodyned signal is proportional to the power in the reference 

15 beam and can be substantially larger than the pressure amplitude of photoacoustic waves 
generated by blood bolus 40 responsive to illumination by light pulse 42. The heterodyned 
signal can be detected by appropriate filtering of signals generated by at least one transducer 
26 using techniques well known in the art. 

The inventors estimate that a reference acoustic beam that focuses about lmW of 

20 power in a blood bolus having a volume of about 1 mm 2 , focuses sufficient energy in the 
blood bolus to couple the photoacoustic response to the reference acoustic beam. The 
coupling of two acoustic signals in the human body to generate a heterodyned signal is 
discussed in a book entitled "Physical Ultrasonics" by Robert T. Beyer and Stephen V. 
Letcher, Academic Press 1969, the disclosure of which is incorporated herein by reference. 

25 Calculations in the book indicate possible gain factors. 

It is to be noted that using a heterodyning detection technique, in accordance with an 
embodiment of the present invention, enables simultaneous excitation and subsequent 
simultaneous detection of photoacoustic waves generated in bolus 40 by light at a plurality of 
different wavelengths. For example, assume that blood bolus 40 is simultaneously illuminated 

30 with light at three different wavelengths in order to determine glucose concentration in the 
blood bolus. Assume further that the pulse rate of the light is different for each wavelength 
and that during illumination a reference acoustic wave is focused on the bolus. The 
photoacoustic response of the blood bolus will result in generation of a strong photoacoustic 
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signal at a different beat frequency for each of the three wavelengths. The signals at each of 
the frequencies can be simultaneously detected and identified using appropriate filtering and 
signal processing techniques known in the art. 

It is to be further recognized that a heterodyning detection technique, in accordance 
5 with an embodiment of the present invention, can be used to determine a location of the 
source of photoacoustic waves generated in a body. Energy in a reference acoustic beam can 
be focused to a relatively small focal volume at a known location. Photoacoustic signals 
generated in the body that are characterized by a frequency equal to a beat frequency must 
have originated in the reference beam focal volume. 

10 Figs. 2A-2C schematically show glucometer 20', similar to glucometer 20, measuring 

concentration of glucose in a person's blood by sensing ultrasound reflections from the blood, 
in accordance with an embodiment of the present invention. 

In Fig. 2A, as in Fig. 1 A, controller 30 controls acoustic transducer 26 to transmit at 
least one ultrasound pulse 32 into body part 24 to locate a bolus of blood 40 in at least one 

15 blood vessel 22 and, optionally to perform acoustic calibration measurements. Controller 30 
uses the intensity of at least one ultrasound pulse 32 and corresponding reflected ultrasound 
pulses 38 from interfaces between walls 41 of blood vessel 22 and blood bolus 40 to 
determine an acoustic energy reflectance for the interfaces of the blood vessel walls and blood 
bolus 40. 

20 In Fig. 2B, controller 30 controls light source 28 to illuminate blood bolus 40 with at 

least one pulse 42 of light having a wavelength useful for assaying glucose in the bolus. 
Energy absorbed from at least one light pulse 42 by blood bolus 40 and blood vessel walls 41 
cause changes in the acoustic impedance of the bolus and the blood vessel walls. The 
impedance changes result in a change of the acoustic energy reflectance at the interfaces 

25 between bolus 40 and blood vessel walls 41. The amount of energy absorbed by bolus 40 
from at least one light pulse 42, and as a result, the change in acoustic energy reflectance of 
bolus 40, is a function of concentrations of glucose and other substances in the blood that 
absorb light at the wavelength of at least one light pulse 42. 

If "R" represents the acoustic energy reflectance of bolus 40, the magnitude of the 

30 change in reflectance, "AR", can often be approximated by AR « K((3AT) 2 . In the expression 
for AR, K is a constant of proportionality and AT is a temperature difference between bolus 40 
and blood vessel walls 41 resulting from absorption of energy from at least one light pulse 42. 
Assuming that the temperature difference AT is generated mostly from change in temperature 
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of blood bolus 40, AR can be written AR « K[B(|3/Hp)(I 0 Te- ad )ab] 2 where B is a constant of 
proportionality. In the expressions for AR, p is the density of blood bolus 42 and the other 
symbols have the same definitions defined above for the expression for P. 

In accordance with an embodiment of the present invention, a change in acoustic 
5 energy reflectance AR of blood bolus 40 resulting from "heating" by at least one light pulse 42 
is measured by reflecting ultrasound waves from bolus 40 after the bolus is illuminated by the 
at least one light pulse. Changes in AR are measured at each of a plurality of different 
wavelengths. At each wavelength an absorption coefficient otb in blood for light at the 
wavelength is determined from the measured change in reflectance. The determined 

10 absorption coefficients are processed similarly to the way in which absorption coefficients 
determined from amplitudes of photoacoustic waves are processed to determine the 
concentration of glucose in the blood. 

In some embodiments of the present invention, controller 30 adjusts values of AR 
using results of the acoustic calibration measurements to remove biases in the values 

15 introduced by acoustic properties of tissue in body part 24 and coupling of at least one 
transducer 26 to the body part. Controller 30 uses the adjusted values of AR to determine 
glucose concentration. 

Fig. 2C schematically shows glucometer 20' measuring acoustic energy reflectance at 
interfaces of blood bolus 40 and blood vessel walls 41 after the bolus is heated by at least one 

20 light pulse 42. In Fig. 2C controller 30 controls at least one ultrasound transducer 26 to 
transmit at least one ultrasound pulse 46 that is incident on bolus 40. Energy from at least one 
ultrasound pulse 46 is reflected from interfaces of bolus 40 and blood vessel walls 41 and 
propagates back to ultrasound transducer 26 as ultrasound pulses 47 of reflected ultrasound 
waves. Transducer 26 generates signals responsive to reflected ultrasound pulse 47, which 

25 signals are transmitted to controller 30. Controller 30 processes the signals to determine an 
acoustic energy reflectance for interfaces of bolus 40 and blood vessel walls 41 using methods 
known in the art. Controller 30 compares the determined reflectance to the reflectance at the 
interfaces prior to illumination by at least one light pulse 42 to determine a change in the 
reflectance. 

30 Preferably, ultrasound pulse 46 is transmitted following a delay from a time at which 

bolus 40 is illuminated by at least one light pulse 42 that is sufficiently long so that at the 
transmission time of ultrasound pulse 46 photoacoustic waves generated by at least one light 
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pulse 42 have subsided. This assures that photoacoustic waves incident on ultrasound 
transducer 26 do not interfere with measurements of reflected ultrasound pulse 42. 

In some embodiments of the present invention, at least one light pulse 42 comprises a 
train of light pulses and glucometer 20' measures a change in reflectance following each light 
5 pulse in the train of light pulses. Optionally, the pulse width and pulse repetition rate are such 
that between light pulses bolus 40 cools down to a temperature substantially equal to the 
ambient temperature. Optionally, energy reflectance of bolus 40 is measured by reflecting an 
ultrasound pulse 46 from the bolus before and after each light pulse in the train of light pulses 
to determine a change in reflectance resulting from illumination by the light pulse. In some 
10 embodiments of the present invention, the pulse width of the light pulses is about 10 
nanoseconds and the pulse repetition rate is about 10 kHz. In some embodiments of the 
present invention, ultrasound pulses 46 are pulses of ultrasound waves at a frequency of at 
least 10 MHz. 

In some embodiments of the present invention changes in acoustic impedance of blood 

15 bolus 40 are measured to determine blood glucose levels while bolus 40 is being heated by at 
least one light pulse 42 or shortly thereafter. To perform the measurements glucometer 20' 
reflects ultrasound from bolus 40 while the bolus is being heated by light pulse 42, or shortly 
thereafter, during a time period in which photoacoustic waves are being generated in bolus 40 
responsive to illumination by the light pulse. 

20 Ultrasound waves reflected from bolus 40 while it is expanding following absorption 

of energy from light pulse 42 often carry a very definite and relatively easily identified 
signature - a very large Doppler shift. During expansion, while a photoacoustic wave is being 
generated in blood bolus 40 a surface of a volume of blood in bolus 40 expands at a speed 
close to the speed of sound. As a result, sound waves reflected from bolus 40 while it is being 

25 heated by light pulse 42, are generally characterized by very large Doppler shifts. 

Figs. 3 A and 3B show a glucometer 21, similar to glucometer 20, operating to assay 
glucose in a person's blood by reflecting ultrasound waves from blood bolus 40 while the 
blood bolus is being heated by light that is absorbed by blood. 

In Fig. 3 A, as in Fig. 2 A, glucometer 21 locates at least one blood vessel 22 and blood 

30 bolus 40 and, optionally, performs acoustic calibration measurements. In Fig. 3B controller 30 
controls light source 28 to illuminate bolus 40 with a light pulse 42 and while bolus is being 
illuminated by the light pulse, controls ultrasound transducer 26 to transmit a pulse 51 of 
ultrasound that is focused on the blood bolus. As a result of illumination of bolus 40 by light 
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pulse 42, bolus 40 expands rapidly. The surface of the bolus moves outward from the bolus at 
speeds close to the speed of sound and generates a photoacoustic wave. Energy from 
ultrasound pulse 51 is reflected from the rapidly moving surface of bolus 40 in an ultrasound 
pulse 54 having a large Doppler shift generated by the velocity at which the surface moves. 
5 The intensity of Doppler shifted pulse 54 is a function of the energy reflectance of the surface 
of bolus 40 which in turn is a function of a difference between the absorption coefficient of 
light in blood bolus 40 and the absorption coefficient of tissue surrounding the bolus. 

In accordance with an embodiment of the present invention, glucometer 21 measures 
intensity of Doppler shifted waves reflected from a blood bolus while the blood bolus is being 

10 illuminated with light for different wavelengths of the illuminating light The measurements 
are used to determine absorption coefficients for blood at the wavelengths and a concentration 
of glucose in the blood. In some embodiments of the present invention, the measurements are 
corrected responsive to acoustic calibration measurements and the corrected measurements are 
used to determine glucose concentration. 

15 Fig. 4 schematically shows a glucometer 60 measuring blood glucose by sensing 

changes in velocity of sound in blood, in accordance with an embodiment of the present 
invention. 

Glucometer 60 is similar to glucometer 20. However, whereas glucometer 20 is shown 
with at least one ultrasound transducer 26, which may be a plurality of ultrasound transducers 

20 located at different positions on body part 24, glucometer 60 is explicitly shown with at least 
one ultrasound transducer 62 and at least one ultrasound transducer 64. At least one 
ultrasound transducer 62 and at least one ultrasound transducer 64 are located substantially 
opposite each other on body part 24. 

As in glucometer 20 and the methods for assaying blood glucose described in the 

25 discussions of Figs. 1 A-3B, glucometer 60 locates at least one blood vessel 22 and a blood 
bolus 40 in the blood vessel using ultrasound and, optionally, performs acoustic calibration 
measurements. In addition glucometer 60 determines an average speed of sound in tissue 
between ultrasound transducers 62 and 64. Different techniques known in the art may be used 
to determine the speed of sound in the tissue. For example glucometer 60 may determine a 

30 transit time between transducers 62 and 64 and use the transit time to determine a speed of 
sound in the tissue. Preferably, glucometer 60 determines an average speed of sound for bolus 
40 responsive to ultrasound reflections from walls 41 of blood vessel 22 using techniques 
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known in the art. Methods for measuring transit times of ultrasound pulses through tissue are 
described in US Patent 4,059,010, the disclosure of which is incorporated herein by reference. 

In Fig. 4 controller 30 controls at least one ultrasound transducer 62 to transmit a 
modulated beam 68 of ultrasound that passes through bolus 40 and is incident on at least one 
5 ultrasound transducer 64. In some embodiments of the present invention, the amplitude of 
beam 68 is modulated. In Fig. 4, by way of example, beam 68 has an amplitude modulation 
envelope 70 shown to the left of the beam. 

At a time to while beam 70 is being transmitted, controller 30 controls light source 28 
to illuminate at least one bolus with a pulse 42 of light having a wavelength that is absorbed 

10 by glucose. Bolus 40 and tissue in body part 24 surrounding bolus 40 absorb energy from 
light pulse 42. The absorbed energy changes the speed of sound in bolus 40 and the 
surrounding tissue. As a result, the transit time of ultrasound from at least one transducer 62 
through bolus 40 to at least one transducer 64 changes. 

It is assumed that a pulse width of light pulse 42 and a rise time during which a 

15 maximum change in the speed of sound in tissue illuminated by the light pulse is established 
are short compared to a transit time of sound through bolus 40. It is further assumed that once 
the maximum change is established, the change is relatively stable during a period of time 
equal to a transit time of sound through tissue in the neighborhood of blood bolus 40 that is 
illuminated by light pulse 42. {i.e. It is assumed that the rise time of the changes in the speed 

20 of sound caused by light pulse 42 is short compared to the transit time of sound through bolus 
40. It is further assumed that the decay time of the changes is long compared to the transit 
time through the volume of tissue illuminated by light pulse 42 in a neighborhood of bolus 
40.) 

Let D be the path length from at least one transducer 62 to at least one transducer 64 
25 and let distance from at least one transducer 62 along the path length be represented by a 
variable x. Let Cq(x) be the speed of sound at position x before illumination with at least one 
pulse 42 and let TTq be a transit time between transducers 62 and 64 before illumination. A 
portion of beam 68 located at a distance x at time to will have a transit time TT(x) that differs 
from . . the transit time TTq . by an amount 

D D D 

30 ATT(x)= J[l/C (x)-l/C(x)]dx= J[(C(x)-C (x))/C (x)C(x)]dx= J[AC(x)/C (x)C(x)]dx . In 

XX X 

the expressions for ATT(x), Cq(x) is the speed of sound at location x before location x is 
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illuminated by light pulse 42 and C(x) is the speed of sound at the location after illumination 
by light pulse 42. 

The derivative of TT(x) with respect to x is equal d[ATT(x)]/dx = AC(x)/Cq(x)C(x). 
AC(x) is generally small compared to Cq(x) and assuming that variations in Cq(x) as x 
5 changes are relatively small, the expression for d[ATT(x)]/dx may be approximated as 
AC(x)/C 2 , where C is an average speed of sound over the path length from at least one 
transducer 62 to at least one transducer 64. It is to be noted that an error in Cq(x) of a 
magnitude on the order of AC(x) will not substantially affect the value for d[ATT(x)]/dx and 
leads to a percent error in d[ATT(x)]/dx on the order of 100[AC(x)/C]%. In situations for 
10 which variations in Cq(x) are not small over the path length, Co(x) can be estimated from 
known speeds of sound in tissue types or determined using ultrasound techniques known in 
the art. 

Values for the derivative d[ATT(x)]/dx and thereby AC(x) at a location x can be 
estimated by using a difference between changes in transit times for two portions of beam 68 
15 that are located in a neighborhood of location x at time tq. In particular d[ATT(x)]/dx and 
thereby AC(x) can be evaluated for locations of x inside bolus 40. 

For example, assume that a minimum 72 in amplitude of beam 68 is located at time to 
at a position xj, which borders bolus 40 on a side of the bolus facing transducer 62. The 
transit time of minimum 72 in comparison to transit time TTq is changed by an amount 

D 

20 ATT(xi) = j[AC Q (x)/C 0 (x)C(x)]dx . Assume that an amplitude maximum 74 at time to is 

x i 

located at position X2 which borders bolus 40 on a side of bolus 40 facing transducer 64. A 
change in transit time for maximum 74 in comparison to transit time TTo is equal to 

D 

ATT(x 2 ) = J[AC 0 (x)/C 0 (x)C(x)]dx . Then [ATT(x 2 ) -ATT(x 1 )]/(x 2 -xi) m d[ATT(x b )]/dx = 

X 2 

AC(xb)/c2, where x b is a location inside bolus 40. It is therefore seen that a value for AC(x) 
25 for blood bolus 40 can be determined using changes in transit times of minimum 72 and 
maximum 74. 

Since the magnitude of AC(xb) is a function of an amount of energy absorbed by bolus 
40 from light pulse 42 the value AC(x b ) can be used to determine the absorption coefficient in 
blood for the wavelength of light in at least one light pulse 42. 
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For example, a change in the speed of sound in bolus 40 caused by energy absorbed by 
the bolus from at least one light pulse 42 is generally approximately proportional to a change 
in temperature in the bolus caused by the absorbed energy. For water, for example, in a 
temperature range from 0° to 100°, AC » 4.9 AT. A change in the speed of sound in bolus 40 
5 as a result of illumination by at least one light pulse 42 may therefore be written (for small 
temperature changes) AC = B[(l/Hp)(I 0 Te~ a d)ctb], where B is a constant of proportionality. 
From the expression for AC it is seen that AC is proportional to the absorption coefficient, a^, 
in blood for light in at least one light pulse 42 and that the absoiption coefficient for the light 
can be determined from AC. 

10 In accordance with an embodiment of the present invention, glucometer 60 determines 

changes in speed of sound in bolus 40, using changes in transit times, for a plurality of 
wavelengths of light. Glucometer 60 uses the determined changes in the speed of sound to 
determine an absorption coefficient for each of the wavelengths and uses the absorption 
coefficients to determine concentration of glucose in the blood. In some embodiments of the 

15 present invention, determined changes in speed of sound in bolus 40 are adjusted responsive 
to acoustic calibration measurements and the adjusted measurements are used to determine 
glucose concentration. 

In some embodiments of the present invention controller 30 controls at least one 
transducer 62 to continuously transmit modulated beam 68 through bolus 40 and to transmit a 

20 plurality of light pulses 42 through bolus 40. Transit times are measured shortly before and 
during each light pulse in the train of light pulses. 

In some embodiments of the present invention, glucometer 60 determines blood 
glucose by sensing Doppler frequency shifts in ultrasonic waves transmitted through blood 
while the blood is continuously exposed to light absorbed by glucose. 

25 Referring again to Fig. 4, as in the case when glucometer 60 operates to assay blood 

glucose from transit times, when operating to assay blood glucose from Doppler shifts in 
transmitted ultrasound, controller 30 controls transducer 62 to transmit a modulated beam 68 
of ultrasound through the bolus. However, whereas when operating in a "transit time" mode, 
at a time to, glucometer 60 illuminates bolus 40 with a light pulse 42 that is relatively short 

30 with respect to a transit time of sound through a region illuminated by the light pulse, when 
operating in a "Doppler shift" mode, glucometer 60 illuminates bolus 40 with a light pulse 42 
that is relatively long with respect to the transit time. As a result, from time to, bolus 40 and a 
region surrounding tissue absorb energy substantially continuously as portions of beam 68 
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transit the bolus and illuminated region. The absorbed energy causes the speed of sound in 
bolus 40 and the region of illuminated surrounding tissue to change continuously during the 
transit of the portions of beam 68. 

As noted above, a sound wave having a frequency f moving through a region in which 
5 the speed of sound is changing, undergoes a Doppler shift "Af" in frequency, which may be 
written Af =f(Ax/C(x) 2 )(dC(x)/dt). In the expression for Af, Ax is the path length of the sound 
wave through the region, x is the location of Ax in the region, dC(x)/dt is the time derivative 
of C(x) at location x. It is assumed that the speed of sound and its derivative are constant over 
the path length Ax. 

10 Therefore, as portions of beam 68 pass through bolus 40 and tissue in the 

neighborhood of bolus 40 illuminated by light pulse 42, the portions undergo a Doppler 
frequency shift. The magnitude of the Doppler shift for a portion of beam 68 received by 
transducer 64, which at time tQ was located at a distance x from transducer 62, is given by 
D 

Af(x) = f J[dC(x)/dt][l/C(x) 2 ]dx . From an analysis similar to that used in the discussion of 
x 

15 determination of glucose levels from changes in transit times, it can be shown that 
d[Af(x)]/dx, and thereby [l/C(x) 2 ][dC(x)/dt)], at a location x can be estimated by using a 
difference between Doppler shifts for two portions of beam 68 that are located in a 
neighborhood of x at time tg. 

Since dC/dt is a function of the rate of absorption of energy from light pulse 42, which 

20 in turn is a function of the absorption coefficient of the light in blood, a measurement of Af 
can be used to determine the absorption coefficient of the light at location x. 

By way of example, for bolus 40 assume that a change in speed of sound, AC, as a 
function of energy absorbed by the bolus is given by the expression for AC above and that C 
and dC/dt are substantially constant in the bolus. Then an expression for dC/dt can be 

25 determined by differentiating the expression for AC with respect to x so that dC/dt = 
B[(l/Hp)(I 0 e" ad )ab] and Af = ffl[(l/c2Hp)(I 0 e-ad) ab ]x, where X is the path length through 
bolus 40 and C is an average value of the speed of sound in the bolus. From the expression for 
Af it is seen that the absorption coefficient ccb can be determined from Af. 

In accordance with an embodiment of the present invention, glucometer 60 measures 

30 Af and determines the absorption coefficient of light in bolus 40 at a plurality of wavelengths, 
and determines the concentration of glucose in blood bolus 40 using the determined 
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absorption coefficients. In some embodiments of the present invention, determined values for 
Af are adjusted responsive to acoustic calibration measurements and the adjusted 
measurements are used to determine absorption. 

In some embodiments of the present invention, a glucometer in accordance with an 
5 embodiment of the present invention assays glucose in a region of tissue by transmitting a 
pulse of ultrasound that passes through the region more than once after the region has been 
illuminated with a pulse of light that is absorbed by glucose. The process is schematically 
illustrated in Figs. 5A and 5B in which a glucometer 80, similar to glucometer 20, is shown 
assaying glucose in a region 90 of tissue in a body part 92. Region 90 is located between bone 

10 tissue 94 and transducer 26 and the assay is performed by transmitting at least one pulse 96 of 
ultrasound that traverses region 90 a plurality of times. The assay is not specific to blood in 
region 90 but is a function of glucose in interstitial fluid and other tissues in the region. As 
with measurements acquired with other glucometers in accordance with embodiments of the 
present invention, in some embodiments of the present invention, measurements acquired 

15 with glucometer 80 that are used to determine glucose concentrations are corrected responsive 
to acoustic calibration measurements. 

In Fig. 5 A controller 30 controls transducer 26 to transmit at least one pulse 96 of 
ultrasound into tissue region 90. Because of differences in acoustic impedance between tissue 
region 90 and bone 94 and between tissue region 90 and transducer 26 and/or air bounding 

20 region 90, ultrasound pulse 96 is repeatedly reflected back and forth between bone 94 and 
transducer 26. As a result, at least one ultrasound pulse 96 repeatedly travels back and forth 
between transducer 26 and bone tissue 94. (In Fig. 5A at least one ultrasound pulse 96 is 
schematically shown traversing tissue region 90 only twice to prevent clutter.) As it travels, at 
least one ultrasound pulse 96 loses energy and is attenuated as a result of absorption in tissue 

25 in region 90 and partial reflection at the interfaces between tissue in region 90 and bone 94 
and tissue in region 90 and transducer 26 or air. 

In response to rebounding ultrasound pulse 96, transducer 26 generates a train of 
signals having a period substantially equal to a round-trip time of ultrasound pulse 96 between 
bone 94 and transducer 26 and decreasing signal amplitude. Controller 30 uses the signals to 

30 determine a transit time "TT(n)" for a given plurality of * V 9 round trips of ultrasound pulse 96 
through tissue region 90. Assuming that the distance between transducer 26 and bone 94 is D, 
then TT(n) = nD/C, where C is the speed of sound in tissue in region 90. Assuming that 
accuracy of measurement of time is constant, as n increases accuracy of measurement of C 



30 



WO 02/15776 PCT/IL01/00740 

increases. In some embodiments of the present invention, the decrease in signal amplitude is 
used to provide an estimate of an absorption coefficient for ultrasound in tissue region 90. 

In Fig. 5B controller 30 controls light source 28 to illuminate tissue region 90 with a 
pulse 42 of light that is absorbed by glucose, after which glucometer 80 measures a transit 
5 time TT*(n) for n round trips of at least one pulse of ultrasound pulse through tissue region 90. 
Energy absorbed from light pulse 42 causes a change AC in the speed of sound in tissue in 
region 90 and TT'(n) = nD/(C+AC). In some embodiments of the present invention controller 
30 determines AC from AC = C{TT(n)/ TT'(n)-l} "and uses AC to assay glucose in tissue 
region 90. In some embodiments of the present invention controller 30 determines a change in 

10 the absorption coefficient for ultrasound that is caused by energy absorbed from light pulse 
42, and uses the change in absorption coefficient to assay glucose in tissue region 90. 

In Fig. 5A and 5B ultrasound pulse 96 is shown being "bounced" between two 
impedance interfaces, a bone 94 and tissue 90 interface and a tissue 90 and ultrasound 
transducer 26 interface. In accordance with embodiments of the present invention, ultrasound 

15 pulses can be bounced back and forth between other types interfaces at which acoustic 
impedance changes to measure characteristics of material between the interfaces. For 
example, in accordance with an embodiment of the present invention, ultrasound pulses are 
bounced back and forth between two ultrasound transducers on opposite sides of a tissue 
region to determine a change in sound velocity and/or absorption coefficient in the tissue 

20 region. In some embodiments of the present invention, ultrasound pulses are bounced back 
and forth between an ultrasound transducer on one side of a tissue region and an interface 
between the tissue region and air on an opposite side of the tissue region to determine a 
change in sound velocity and/or absorption coefficient in the tissue region. In some 
embodiments of the present invention, ultrasound pulses are bounced between blood vessel 

25 walls to determine a change in sound velocity and/or absorption coefficient in blood in the 
blood vessel. 

Other variations of the method, in accordance with embodiments of the present 
invention, in which ultrasound is bounced between acoustic impedance interfaces to measure 
characteristics of material between the interfaces will readily occur to persons of the art. 
30 For example, in some embodiments of the present invention, speed of sound in blood 

in a blood vessel is determined by transmitting a first pulse of ultrasound from a transducer so 
that it is reflected back to the transducer from a "near" wall of a blood vessel facing the 
transducer. A second ultrasound pulse is transmitted from the transducer so that it is reflected 
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back to the transducer from a far wall of the blood vessel. The first pulse does not travel 
through blood in the blood vessel to travel back and forth to the transducer. The second pulse 
on the other hand, in travelling back and forth to the transducer travels through blood in the 
blood vessel over a path length substantially equal to twice the diameter of the blood vessel. 
5 In accordance with an embodiment of the present invention, the diameter of the blood vessel 
is determined from ultrasound measurements of the location and size of the blood vessel. The 
speed of sound in blood in the blood vessel can then be determined from the back and forth 
transit times to the transducer of the first and second ultrasound pulses. A change in the speed 
of sound in the blood resulting from illumination of the blood is determined by measuring 

10 back and forth transit times from near and far walls of the blood vessel before and after 
illumination. The change in the speed of sound determined in this manner is used, as in other 
embodiments of the present invention, to assay blood glucose. 

In some embodiments of the present invention, glucometer 80 determines glucose 
concentration in tissue region 90 by determining a "resonant" frequency for ultrasound 

15 transmitted into region 90 for which ultrasound reflected from bone 94 has either a minimum 
or a maximum. 

Resonant frequencies are defined by a matching condition between the wavelength of 
the ultrasound and the thickness of region 90 (i.e. the distance between transducer 26 and 
bone 94) is met. For example, reflection maxima and minima (and respectively transmission 

20 minima and maxima) occur for an incidence angle 0 of ultrasound on bone 90 at frequencies 
"f 1 for which the ultrasound in tissue region 90 has a wavelength X = C/f =* 4Dcos8/n. In the 
expression for X, n is equal to an odd and even integer respectively for reflection maxima and 
minima. From the expression it is seen that values for the resonant frequencies are dependent 
upon the speed of sound C in tissue region 90. 

25 In accordance with an embodiment of the present invention, glucometer 80 determines 

a resonant frequency for sound transmitted into region 90 and then illuminates the region with 
light absorbed by glucose. Energy absorbed from the light causes a change in the speed of 
sound in region 90 and results in a shift in the resonant frequency. The magnitude of the shift 
is measured and used to determine a magnitude for the change in the speed of sound, which 

30 magnitude is used to determine an absorption coefficient for the light in tissue in region 90. It 
is to be noted that whereas embodiments of the present invention have been described for 
assaying glucose, methods and apparatus, in accordance with embodiments of the present 
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invention, are useable to assay substances other than glucose and to assay substance in tissues 
other than blood. 

For example, whereas the method described with reference to Figs. 5A and 5B, in 
which ultrasound is passed through a tissue region a plurality of times, is used to assay 
5 glucose in the tissue, the method can be used to determine concentrations of components other 
than glucose in tissue. For example, at a bone-soft tissue interface, acoustic energy reflectance 
is relatively high. Ultrasound waves can be transmitted into bone tissue and repeatedly 
bounced back and forth between bone-soft tissue interfaces to determine the speed of sound 
and acoustic absorption coefficient in bone. The bone can then be "excited" with 
1 0 electromagnetic radiation, such as optical, radio frequency or microwave radiation, which is 
absorbed by bone calcium. A change in the speed of sound and/or acoustic absorption 
coefficient resulting from energy absorbed from the exciting radiation is then measured. The 
measured change, in accordance with an embodiment of the present invention, can be used to 
determine bone density. 

15 By way of another example of assaying substances other than glucose, in accordance 

with embodiments of the present invention, by choosing an appropriate set of wavelengths for 
light used to illuminate regions of a person's body, concentrations of cholesterol in the 
regions can be determined. Light at wavelengths between 1.70-1.80 microns is absorbed by 
cholesterol and light at these wavelengths can be used in the practice of the present invention 

20 to determine cholesterol concentrations in regions of a person's body. 

Methods of detennining concentration of a substance in a body, in accordance with an 
embodiment of the present invention provide the concentration of the substances as a function 
of position in the body. In some embodiments of the present invention a spatial map of 
concentrations of a substance in a body is used to image features and/or components of the 

25 body. 

Figs. 6A and 6B schematically show a "glucometer" 100, hereinafter referred to as 
"assay imager 100", used, by way of example, to assay and image cholesterol in a person's 
body in accordance with an embodiment of the present invention. Assay imager 100 
comprises a visual display 102. In Fig. 6A, controller 30 controls at least one transducer 26 to 
30 locate a region 36 of at least one blood vessel 22. In Fig. 6B assay imager 100 then 
illuminates region 36 with light at wavelengths suitable for assaying cholesterol and 
determines concentrations of cholesterol in region 36 using methods similar to those described 
above for determining concentrations of glucose. Controller 30 then controls transducer 26 to 
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locate other regions of at least one blood vessel 22 and proceeds to determine concentrations 
of cholesterol in these other regions. A map 110 of cholesterol concentration as a function of 
position for regions scanned by assay imager 100 is displayed on visual display 102. At least 
one blood vessel 22 is compromised by plaque deposits 104 and 106 at two locations. The 
5 locations of plaque deposits 104 and 106 are visible on cholesterol concentration map 110 
shown on display 1 02. 

By way of another example, assay imager 100 may also be used to image tumors. 
Tumorous growths generally have an unusually high concentration of blood vessels in the 
growths and in tissue surrounding the growths. By imaging a region of tissue with assay 

10 imager 100 using light at a wavelength that is strongly absorbed by blood, in accordance with 
an embodiment of the present invention, tissue having an unusually high density of blood 
vessels can be clearly contrasted against tissue having normal blood vessel density. As a 
result, tumor tissue in the imaged region can often be clearly contrasted with healthy tissue. 

In some embodiments of the present invention, an "assay image" of a region of the 

15 body responsive to concentration of a substance in the region that is provided by an assay 
imager is overlaid on an image of the region provided by another imaging modality, such as 
for example an MM or CT image. In some embodiments of the present invention, the image 
on which the assay image is overlaid is an ultrasound image. Optionally the ultrasound image 
is acquired using the same ultrasound detectors that are used to acquire the assay image. By 

20 using the same ultrasound detectors for acquiring both the assay image and the ultrasound 
image on which the assay image is overlaid the assay image is automatically registered to the 
ultrasound image. 

In the description and claims of the present application, each of the verbs, "comprise" 
"include" and "have", and conjugates thereof, are used to indicate that the object or objects of 

25 the verb are not necessarily a complete listing of members, components, elements or parts of 
the subject or subjects of the verb. 

The present invention has been described using detailed descriptions of embodiments 
thereof that are provided by way of example and are not intended to limit the scope of the 
invention. The described embodiments comprise different features, not all of which are 

30 required in all embodiments of the invention. Some embodiments of the present invention 
utilize only some of the features or possible combinations of the features. Variations of 
embodiments of the present invention that are described and embodiments of the present 
invention comprising different combinations of features noted in the described embodiments 
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will occur to persons of the art. The scope of the invention is limited only by the following 
claims. 



35 



4 



t 



WO 02/15776 PCT/IL01/00740 

CLAIMS 

1. A method for assaying a component of a localized region of interest in a body 
comprising: 

5 a) illuminating the region with at least one pulse of radiation having a wavelength at 

which the radiation is absorbed by the component to generate a change in an acoustic property 
of the region; 

b) transmitting ultrasound so that it is incident on the region; 

c) measuring at least one effect of the change on the incident ultrasound; 

10 d) using the measured at least one effect to determine an absorption coefficient for the 

radiation in the region; and 

e) using the determined absorption coefficient to determine concentration of the 
component in the region. 

15 2. A method according to claim 1 and comprising illuminating the region with at least 
one pulse of radiation at at least one other wavelength, repeating b-d and using the determined 
absorption coefficient for the at least one other wavelength to determine concentration of the 
component in the region. 

20 3. A method according to claim 1 or claim 2 wherein an effect of the at least one effect 
comprises a change in acoustic energy reflectance of the region. 

4. A method according to claim 3 wherein measuring the change in reflectance 
comprises: 

25 reflecting at least one first pulse of ultrasound from the region before illumination of 

the region with the at least one radiation pulse; 

determining a first acoustic energy reflectance from the region using the amplitude of 
the at least one reflected pulse; 

reflecting at least one second pulse of ultrasound from the region after illumination of 
30 the region with the at least one radiation pulse; 

determining a second acoustic energy reflectance from the region using the amplitude 
of the second at least one reflected pulse; 

determining a difference between the first and second acoustic energy reflectances. 
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5. A method according to claim 3 wherein measuring the change in reflectance 
comprises: 

measuring acoustic energy reflectance from the region before illumination with the at 
least one radiation pulse; 

5 measuring amplitude of Doppler shifted ultrasound waves reflected from the incident 

ultrasound during illumination of the region with the at least one light pulse; 

using the amplitude of the Doppler shifted ultrasound waves to determine acoustic 
energy reflectance that characterizes the region during illumination with the at least one 
radiation pulse; and 

10 determining a difference between the reflectance before illumination and during 

illumination. 

6. A method according to any of claims 1-5 wherein an effect of the at least one effect 
comprises a change in speed of sound in the region. 

15 

7. A method according to claim 6 wherein measuring the change in speed of sound 
comprises: 

transmitting an ultrasound wave from a first transducer towards a second transducer so 

i 

that the ultrasound wave traverses the region; 
20 illuminating the region while the ultrasound wave is present in the region; 

determining a first transit time from the first to the second transducer of a first portion 
of the ultrasound wave, which first portion is present at a first location in the region while the 
region is illuminated; 

determining a second transit time from the first to the second transducer of a second 
25 portion of the ultrasound wave which second portion is present at a second location in the 
region at a time at which the first portion is located at the first location; and 

using the difference in the transit times to determine a change in the speed of sound 
through the region at a location in the region between the first and second locations. 

■ 

30 8. A method according to claim 6 wherein the region is sandwiched between a first and a 
second acoustic interface at which ultrasound is partially reflected and wherein measuring the 
change in the speed of sound comprises: 
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measuring before illumination a difference between a transit time of ultrasound 
between a first and second location that is reflected from the first interface and a transit time 
of ultrasound energy between the first and second locations that is reflected from the second 
interface; 

5 determining a speed of sound before illumination for the region from the transit time 

difference; 

measuring after illumination a difference between a transit time of ultrasound between 
a first and second location that is reflected from the first interface and a transit time of 
ultrasound energy between the first and second locations that is reflected from the second 
10 interface; 

determining a speed of sound after illumination for the region from the transit time 
difference; and 

determining the change in the speed of sound from the determined speeds of sound 
before and after illumination. 

15 

9. A method according to claim 6 wherein the region is sandwiched between two 
acoustic interfaces at which ultrasound is partially reflected and wherein measuring the 
change in the speed of sound comprises: 

a) introducing a pulse of ultrasound into the region before illumination so that the 
20 pulse bounces back and forth between the interfaces and repeatedly traverses the region; 

b) sensing energy from the pulse that exits the region through one of the interfaces for 
each of a plurality of back and forth txaversals of the pulse through the region; 

c) determining a speed of sound in the region from the times at which the energy is 

sensed; 

25 d) introducing a pulse of ultrasound into the region after illumination so that the pulse 

bounces back and forth between the interfaces and repeatedly traverses the region and 
repeating steps b and c; and 

e) determining a difference between the determined speeds of sound. 

■ 

30 10. A method according to claim 6 wherein the region is sandwiched between two 
acoustic interfaces at which ultrasound is partially reflected and wherein measuring the 
change in speed of sound comprises: 
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transmitting ultrasound which is incident on the region prior to illumination of the 
region with the at least one radiation pulse; 

determining a first frequency for which transmission of the ultrasound through the 

region is . a maximum or a minimum; 
5 transmitting ultrasound which is incident on the region after illumination of the region 

with the at least one radiation pulse; 

determining a second frequency for which transmission of ultrasound through the 
region is a maximum or a minimum and for which the wavelength of propagation of 
ultrasound in the region for the first and second frequencies is substantially the same; and 
10 using a difference between the first and second frequencies to determine the change in 

the speed of sound in the region. 

11. A method according to claim 6 wherein the region is sandwiched between two 
acoustic interfaces at which ultrasound is partially reflected and wherein measuring the 

1 5 change in speed of sound comprises: 

transmitting ultrasound which is incident on the region prior to illumination of the 
region with the at least one radiation pulse; 

determining a first frequency for which reflection of the ultrasound through the region 

is a maximum or a minimum; 
20 transmitting ultrasound which is incident on the region after illumination of the region 

with the at least one radiation pulse; 

determining a second frequency for which reflection of ultrasound through the region 
is a maximum or a minimum and for which the wavelength of propagation of ultrasound in 
the region for the first and second frequencies is substantially the same; and 
25 using a difference between the first and second frequencies to determine the change in 

the speed of sound in the region. 

12. A method according to any of claims 1-11 wherein an effect of the at least one effect 
comprises a change in the frequency of ultrasound that traverses the region. 

30 

13. A method according to claim 12 wherein measuring the change in frequency 
comprises: 
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transmitting an ultrasound wave from a first transducer towards a second transducer so 
that the ultrasound wave traverses the region; 

illuminating the region with the at least one light pulse while the ultrasound wave is 
present in the region and wherein the at least one light pulse has a pulse width substantially 
5 longer than the transit time of ultrasound through the region; 

determining a first frequency shift in the frequency of a first portion of the ultrasound 
wave, which first portion is present at a first location in the region while the region is 

♦ 

illuminated; 

determining a second frequency shift in the frequency of a second portion of the 
10 ultrasound wave which second portion is present at a second location in the region at a time at 
which the first portion is located at the first location; and 

using a difference in the first and second frequency shifts to determine a frequency 
shift that occurs in the ultrasound at a location in the region between the first and second 
locations. 

15 

14. A method according to any of claims 1-13 wherein an effect of the at least one effect 
comprises a change in absorption of ultrasound along a path length in the region. 

15. A method according to claim 14 wherein the region is sandwiched between two 
20 acoustic interfaces at which ultrasound is partially reflected and wherein measuring the 

change in absorption over the path length comprises: 

a) introducing a pulse of ultrasound into the region before illumination so that the 
pulse bounces back and forth between the interfaces and repeatedly traverses the region; 

b) sensing an amount of energy from the pulse of ultrasound that exits the region 
25 through one of the interfaces for each of a plurality of times that the pulse travels back and 

forth through the region; 

c) determining amounts of energy absorbed from the pulse for each back and forth 

traversal through the region; 

d) introducing a pulse of ultrasound into the region after illumination so that the pulse 
30 bounces back and forth between the interfaces and repeatedly traverses the region and 

repeating b and c; and 

e) determining the change in absorption over the path length using the determined 
amounts of absorbed energy. 
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16. A method according to any of the preceding claims and including determining a 
location of the region. 

5 17. A method according to claim 16 wherein detennining a location for the region 
comprises transmitting ultrasound into the body. 

18. A method according to any of the preceding claims wherein transmitting ultrasound 
incident on the region and measuring an effect of the change on the incident ultrasound is 

10 mediated by at least one ultrasound transducer coupled to the body and wherein the method 
comprises measuring acoustic properties of the coupling between the at least one transducer 
and the body. 

19. A method according to any of the preceding claims wherein the at least one pulse of 
15 radiation comprises a microwave pulse. 

20. A method according to any of the preceding claims wherein the at least one pulse of 
radiation comprises a pulse of RF energy. 

20 21. A method according to any of the preceding claims wherein the at least one pulse of 
radiation comprises a pulse of light. 

22. A method for assaying a component of a localized region of interest in a body 
comprising: 

25 a) illuminating the region with at least one pulse of light having a wavelength at which 

light is absorbed by the component; 

b) sensing photoacoustic waves generated responsive to the light pulse and 
determining locations of their origins; 

c) using intensity of those photoacoustic waves having an origin in the region to 
30 determine an absorption coefficient for the light in the region; and 

d) using the determined absorption coefficient to determine concentration of the 
component in the region. 
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23. A method according to claim 22 wherein the at least one light pulse comprises at least 
one train of light pulses radiated at a pulse repetition rate. 

24. A method according to claim 23 wherein sensing and deteraiining locations for the 
5 photoacoustic waves comprises: 

focussing an acoustic reference beam on the region, which reference beam has a 
frequency that is shifted from the pulse repetition frequency of the at least one light pulse train 
by an offset frequency and an intensity that causes the region to respond non-linearly to 
acoustic stimuli; and 
10 detecting acoustic waves at the offset frequency. 

25. A method according to any of claims 22-24 and comprising illuminating the region 
with at least one pulse of light at at least one other wavelength, repeating b and c and using 
the determined absorption coefficient for the at least one other wavelength to determine 

1 5 concentration of the component in the region. 

26. A method according to claim 24 and comprising illuminating the region with at least 
one other light pulse train characterized by a different wavelength and pulse repetition rate, 
repeating b and c, and using the determined absorption coefficient for the at least one other 

20 light pulse train to determine concentration of the component in the region. 

27. A method according to claim 26 wherein the light pulse train and at least one other 
light pulse train illuminate the region substantially simultaneously. 

25 28. A method according to any of claims 22-27 wherein sensing photoacoustic waves is 
mediated by at least one ultrasound transducer coupled to the body and wherein the method 
comprises measuring acoustic properties of the coupling between the at least one transducer 
and the body. 

30 29. A method according to any of claims 22-28 and including determining a location of 
the region. 
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30. A method according to claim 29 wherein determining a location for the region 
comprises transmitting ultrasound into the body. 



31. A method according to claim 17 or claim 30 wherein the region is located to within an 
5 uncertainty of less than 5 wavelengths of the ultrasound in at least one direction. 

32. A method according to claim 17 or claim 30 wherein the region is located to within an 
uncertainty of less than 3 wavelengths of the ultrasound in at least one direction. 

10 33. A method according to claim 17 or claim 30 wherein the region is located to within an 
uncertainty of about a wavelength of the ultrasound in at least one direction. 

34. A method according to any of claims 21-33 wherein the at least one light pulse 
comprises a pulse of IR light. 

15 

35. A method according to any of claims 21-34 wherein the at least one light pulse 
comprises a pulse of visible light. 

36. A method according to any of the preceding claims and comprising measuring 

■ 

20 acoustic properties of material in the body. 

37. A method according to claim 36 wherein a determined absorption coefficient is 
adjusted responsive to the measured acoustic properties. 

25 38. A method according to any of the preceding claims wherein the body is a living body. 

39. A method according to claim 38 wherein the living body is a human body. 

40. A method, according to claim 38 or claim 39 wherein the region is a bolus of blood in a 
30 blood vessel of the body. 

41 . A method according to claim 40 wherein the component is glucose. 
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42. A method of imaging internal features of a body comprising: 

assaying a component of the features for a plurality of regions in the body according to 
any of claims 1-37; and 

displaying a result of measurements performed in assaying the component as a 
5 function of location of the regions to provide an image of the features. 

43. A method according to claim 42 and comprising overlaying the image provided from 
assaying the component on another image of the region acquired by a different imaging 
modality. 

10 

44. A method according to claim 43 wherein the different imaging modality comprises 
ultrasound imaging. 

45. A method according to claim 44 and comprising using a same at least one ultrasound 
15 transmitter or ultrasound sensor to acquire both the image provided from assaying the 

component and the ultrasound image. 

46. A method according to any of claims 42-45 wherein the body is a living body. 

20 47. A method according to claim 46 wherein the living body is a human body. 

48. A method according to claim 46 or claim 47 wherein the feature is a plaque deposit in 
a blood vessel. 

25 49. A method according to claim 40 or claim 48 wherein the component is LDL 
cholesterol or oxidized LDL cholesterol. 

« 

50. A method according to claim 46 or claim 47 wherein the feature is a tumor. 

30 51 . A method for detecting photoacoustic waves generated in a region of a body by light 
comprising: 

illuminating the region with pulses of the light at a pulse repetition frequency to 
generate photoacoustic waves in the region; 
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heterodyning the generated photoacoustic waves with an acoustic reference beam 
having a frequency that is shifted from the pulse repetition frequency by an offset frequency; 
and 

detecting acoustic waves at the offset frequency. 

52. A method according to claim 51 wherein heterodyning comprises focusing the 
acoustic reference beam in the region so that acoustic energy in the region from the reference 
beam is sufficient to cause the region to respond non-linearly to acoustic stimuli. 
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